Competitive localization of vortex lines and interacting bosons 
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We present a theory for the locahzation of three-dimensional vortex lines or two-dimensional 
bosons with short-ranged repulsive interaction which are competing for a single columnar defect or 
potential well. For two vortices we use a necklace model approach to find a new kind of delocalization 
transition between two different states with a single bound particle. This exchange-delocaUzation 
transition is characterized by the onset of vortex exchange on the defect for sufficiently weak vortex- 
vortex repulsion or sufficiently weak binding energy corresponding to high temperature. We calculate 
the transition point and order of the exchange-delocalization transition. A generalization of this 
transition to arbitrary vortex number is proposed. 
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Introduction. Melting of the Bose glass, the low- 
temperature vortex phase in type-II superconductors 
with columnar defects, remains a subject of constant in- 
terest 0, 01 • The interest is motivated not only by the 
appeal and importance of understanding a basic phase 
transition of the vortex system, but, since the 3D vor- 
tex array is equivalent to a quantum 2D Bose system , 
vortex melting also offers a unique experimentally acces- 
sible model to explore the interplay between disorder and 
interactions in the delocalization transition of the corre- 
sponding strongly correlated 2D quantum system. 

Numerous experimental observations on Bose glass 
melting (see, for example, Refs. revealed a charac- 

teristic kink in the low-field segment of the melting line 
suggesting a change of the melting mechanism. A recent 
experimental study of BSCCO samples with a very low 
density of columnar pins ^] allowed one to identify the 
low-field part of the melting curve as depinning transition 
from a single columnar defect driven by vortex- vortex in- 
teractions. A theoretical study of the interacting boson 
system with a low density of strong defects 0| demon- 
strated the possibility of an intermediate superfluid state 
where condensate and localized bosons coexist. Further- 
more, it was found in Ref. ,5j that interactions suppress 
localization and the increase of the boson density results 
in a sharp delocalization crossover into a state where all 
bosons are delocalized. 

The model that is conventionally used in studies of 
quantum localization can be viewed as an ensemble of (in- 
teracting) particles immersed in a random field that can, 
in principle, localize or accomodate all particles; in other 
words, there is a dense array of pinning centers strug- 
gling to localize dilute, but interacting, particles. Ref- 
erence proposed to take an alternative approach and 
consider quantum particles or vortex lines of high density 
competing for dilute traps or pinning sites. In this Let- 
ter, we extend this approach and explore the regime of 
low particle densities, i.e., a finite number of interacting 



quantum particles. A related quantum system that has 
been studied previously are two electrons interacting in 
the region of a short-range attractive potential 

In this Letter, we investigate the formation of bound 
states in an ensemble of N strongly repulsive particles 
competing for a single attractive potential well. In the 
related vortex system, this corresponds to N vortices 
competing for a single columnar defect, i.e., to a physi- 
cal situation where vortices outnumber columnar defects 
(magnetic fields B = NBg, well exceed the matching field 
B<i, = $o/ai) where a^i is the average defect spacing and 
$0 the flux quantum). We focus on the situation where 
mutual repulsion is strong enough to suppress double- 
occupancy of the potential well and consider transversal 
dimensionalities d < 2, for which a bound state for a sin- 
gle particle in a symmetric potential well always exists 
(in the following, we use the language of either particles 
or vortex lines at our will). The main finding of this 
Letter is a new kind of delocalization phase transition 
driven by the exchange of the single bound particle with 
the — 1 unbound ones. We first derive this result for 
N — 2 and then propose a generalization for arbitrary N. 
Finally we obtain the transition temperature Tde for the 
exchange-delocalization transition in the vortex system 
and discuss the resulting phase diagram Fig. |5| Contrary 
to Ref. where the boson (vortex) density was finite, 
we deal in this work with the genuine thermodynamic 
limit of infinite system size but finite particle (vortex) 
number N corresponding to B, B$ « with N = B/B<j, 
finite. Thus the exchange-delocalization phase transition 
emerges in the limit B ^ and replaces the crossover 
that was found in Ref. 0] for macroscopic vortex density 
B>0. 

The model. We describe a single vortex line in a sam- 
ple of thickness L interacting with an attractive columnar 
defect by the Hamiltonian 

-Hiiriz)] = ^"^dz |ie,(9.r)2 + y^(r(z))| , (I) 
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where (z,r(z)) is the vortex trajectory in 1 + d dimen- 
sions, the magnetic field aUgned with the z axis, and ei 
is the stiffness or tih modulus of a single vortex line; 
in an anisotropic superconductor £/ k, e^gg In k where 
£o = (*i'o/47rA)^ is the characteristic vortex line energy, 
A is the magnetic penetration depth, k = A/^ ^ 1, ^ is 
the coherence length, and e is the anisotropy parameter. 
Vd (r) is the pinning potential from a single columnar de- 
fect positioned at r = 0. Vd{t) falls off exponentially for 
7" > A such that the large scale behaviour of pinned 
vortex lines is well described using a cylindrical pinning 
potential well 

V^(v)^Uj^ for r<hD., Vi)(r) = for r > hu (2) 

with a potential depth C/^j = — eo?'Li/4C'^ and an effective 
radius hu = \J2^ + where ru is the radius of the 
columnar defect 3]. 

In d < 2 dimensions the line is always bound to the 
defect as can be seen via mapping onto the ground-state 
problem of a quantum particle in a potential Voiv) for 
large L Q . Choosing the energy scale such that the un- 
pinned vortex line has a free energy Eq = 0, the free 
energy per length _Ei < of the bound vortex line is ob- 
tained as the ground state energy Ei of the Schrodinger 
equation 

[(TV2e/)V2 - Vnir)] ^Pir) - -i^i^W, (3) 
^1 - C/D/(r/T*), where T* = bDV£i\UD\ (4) 

is a characteristic crossover energy, and f{x) a scaling 
function. For d = 2 and the pinning potential , it has 
the asymptotic behavior f{x) « 1 — C'(x^) for x ^ 1 and 
fix) « x^ exp(-2x2)/2 for a; > 1 3]. 

For A'^ > 1 we add the repulsive vortex interactions 



dzV{\r,{z)-r,{z)\), (5) 
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where V{r) — 2eoKo{r / X), and Kq is the Bessel function. 
Double-occupancy of the defect over an extended length 
is energetically disfavored if Eq ^ which is always 

the case at high enough temperatures. We also focus on 
the regime of large vortex spacing a S> A. 

Exchange- delocalization transition for N—2. We in- 
vestigate the localization behavior for the case N = 2 
making use of a necklace model approach As the 

double-occupancy of the defect is energetically unfavor- 
able, the accessible configurations of the vortex lines con- 
sist of a necklacelike succession of two possible configu- 
rations A and B, see Fig. ^ In the configuration A, 
the line 1 is bound to the defect with the binding free 
energy Ei < and the transversal localization length 
^± — T/ y^\Ei \ei, whereas the line 2 is essentially in the 
unbound state with the free energy Eq = and experi- 
ences rare collisions with the bound line 1. As we assume 



A <C a, the unbound line is exploring the region r > A 
of exponentially weak repulsion, whereas collisions oc- 
cur at r < A. The configuration A ends in an exchange 
event where the endpoints of the unbound line 2 attach 
to the defect again; see Fig. ^ At these exchange points 
the configuration A can connect to the configuration B 
where the roles of the particles are exchanged. 

First, we estimate the energy cost of the localized col- 
lision and exchange events. In the presence of the repul- 
sion lO, each return of line 2 to the bound line 1 will cost 
an additional collision repulsion energy Er which is de- 
termined by optimizing the sum of elastic and repulsive 
energies for a contact of length £r over which the typical 
line spacing is of the order |Ar| ~ A, iJ^ — eiX'^/ir + ^r^o'- 



Ir ~ Xy/ei/so, 



(6) 



and Vr = exp{—Er/T) < 1 defines the Boltzmann-factor 
associated with each collision. Similarly, we estimate the 
energy cost E^x of a localized exchange by optimizing the 
sum of the elastic energy and the loss of binding energy 
\Ei\iex for a contact of length £ex over which the typical 
line spacing is of the order |Ar| ~ A, E^x — eA^/^er + 
£c\Ei\. This gives 



£ex ^ X^ei/\Ei\, Eex ^ X^/I^i 



(7) 



and Vex = exp(— i^ea/T) < 1 is the Boltzmann-factor 
associated with each localized exchange. 

Now we address the statistical mechanics problem of 
summing over all vortex line configurations. Adopting 
a coarse-grained description focusing on scales r > X 
for transversal vortex fluctuations, we discretize the vor- 
tex system into segments of length I — X^£i/T in the 
z-direction. In what follows, we calculate the grand- 
canonical partition sum G{z) — Z{n)z^ where Z{n) 
is the partition sum for a system of length L = nl and z is 
the fugacity. G{E) = G(exp (lE/T)) is the Green's func- 
tion at energy E for the corresponding two-particle quan- 
tum problem. The free energy density / of the system 
is determined by the real singularity zq of G(z) closest 
to the origin by the relation zq = exp{lf /T). If Ga{z) 
and Gb{z) are the partition sums for configurations A 



exchange 




FIG. 1: Two particles binding alternately to a single columnar 
defect. The particle binding to the defect can be exchanged 
in localized events and there are rare collisions between the 
free and the bound particle. 
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and B, respectively, the full partition sum is obtained 
by summing over all alternating configurations Ga, Gb, 
GaVsxGb, GbVsxGa, GAVexGBVexGA,- ■ ■ , Separated by 



particle exchanges with Boltzmann-factor Vg 
up the resulting geometric series we obtain 



Summing 



G(z) = 



Ga + Gb + 2vexGAGB 



1 — vI^GaGi 



'2Ga{z) 



1 - VexGA{zy 



(8) 



where we used Ga = Gb because both configurations 
are related by a mere particle exchange. [Boundary ef- 
fects are irrelevant; we allow either ^ or i? at the ends 
of the defect in ((SJ]. According to Eq. the singu- 
larity determining the free energy of the system is given 
either by the singularity of Ga{z) corresponding to the 
state where the same line is always bound or by the solu- 
tion of 1 = VexGA{z) corresponding to alternating bound 
particles. The exchange- delocalization transitionhetween 
these two states occurs if both singularities occur at the 
same value of z. 

To move further, we calculate the grand-canonical par- 
tition sum Ga^z). In the absence of the interline re- 
pulsion, the canonical partition sum in configuration 
j4, ZA{n) — Zi{n)Z2{n), is a product of the partition 
sum of the bound line 1, Zi{n) = exp{—nlEi/T) [see 
Eq. Q], and of the free line 2, which is attached with 
its end-points to the defect. This restriction leads to 
Z2{n) — p„ exp(— n/i?o/r) where p„ is the probability 
for the return of the unbound line to the defect and 
i?o = is its free energy per unit length. In the ab- 
sence of the repulsion this return probability is given by 
the return probability of a random walk, p„ ~ n"''/^, 
for large n. Then the partition function Ga^z) is re- 
lated to the generating function P(z) = 'YlmVnZ^ for 
these return probabilities by Ga{z) = P(wz) where 
w = eyi^{-lEi/T) = eTcp{-XeiEi/T^). For 1 - z < 1, 

P{z) - (1 - for d < 2 and P{z) ln(l - z) 

for d = 2. Including the Boltzmann-factor Vr in the ran- 
dom walk of the unbound line for each collision with the 
repulsive bound line localized at the defect [sec Eq. jnj] 
leads to a modified generating function Pr{z): 

Pr{z) = P{z) + P{z)VrPr{z) - P{z)Pr{z). (9) 

In this relation the contributions from the repulsion-free 
walks with Boltzmann-factor 1 are subtracted and the 
corresponding term with modified Boltzmann-factor Vr is 
added on the right hand side recursively for each collision. 
With this modification due to collisions, we finally obtain 



P{wz) 



GA{z)^Pr{wz)^——^^^-—- 



for the grand-canonical partition function. Ga{z) has a 
singularity at z — l/w corresponding to a free energy 
per length f = Ei identical to that of a single bound 
particle because the second unbound particle has Eq = Q. 



The function P{wz) diverges upon approaching z = \/w 
and, thus, we find Ga(1/w) = ^r(l) = (1 — Vr)~^ at the 
singularity. 

Now we turn to the exchange-delocalization transition 
determined by the singularities of ©. We have already 
found that the singularity of Ga{z) (at z — l/w corre- 
sponding to the free energy f = Ei) describes, indeed, a 
single bound particle, i.e., a state with always the same 
line bound. We have argued above that there can be a 
real singularity closer to the real axis which is given by 
the solution of the equation 1 = VexGA^z) and which 
corresponds to the phase with exchanging bound parti- 
cle. From the functional form (|10|l of Ga{z), one readily 
verifies that this singularity is indeed the one closer to the 
origin and therefore representing the thermodynamically 
stable phase provided 



> Vex,c = G/{1/w) 1 



(11) 



According to the exchange-delocalization phase 

transition occurs at the critical temperature Tde that is 
obtained from exp(— i^r/Tde) « E^x/Tde where we as- 
sumed that Eex ^ Tde because Tde > T*. Using esti- 
mates © and 13), the asymptotics of the function f{x) 
for d = 2, and X/bn = k/\/2, we arrive at 



T, 



(12) 



for the delocalization transition temperature Tde in the 
vortex system. The transition takes place in the regime 
T > T* where a single vortex is only weakly bound to 
the defect 0; see Fig. |21 Note that both Eqs. Q and 
(|12|l are self-consistent equations for T* and Tde, respec- 
tively, due to the temperature dependence of ^ and A. 
Furthermore, it can be shown by expanding the equa- 
tion 1 = VexGAiz) about the transition point, that the 
exchange-delocalization transition is continuous for all 
d < 2 and of infinite order for d = 2. 



Vortex liquid 




FIG. 2: Schematic phase diagram in the B-T plane. Bm is the 
pristine melting line, Bbg the Bose glass melting line and Bde 
the exchange-delocalization crossover line which terminates 
at S = in a genuine phase transition at temperature Tde 
(circle) . 
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We expect our results to apply to all short-ranged po- 
tential wells and particle interactions that decay faster 
than 1 jr^ for large separations r . More realistic pin- 
ning potentials contain an intermediate 1/r^-bchaviour 
on scales f ^ r ^ A 0, , which slightly changes the 
function /(a;) and thus the exact value of T^e but not the 
universal properties of the delocalization transition. 

General N. We start from the exchange-delocalized 
state where all N vortices share the defect and consider 
the instability with respect to the exclusion of one of the 
vortices from the exchange. To this end we introduce N 
states (analogously to the states A and B for = 2), 
where one vortex is unbound, i.e., excluded from the ex- 
change, whereas the other — 1 vortices share the de- 
fect. Then the necklace is a succession of possible states 
i ~ l,...,N each of which has a grand-canonical par- 
tition sum G^{z) = P{ WNZ)/ {1 Vr{N)]P(wNz)) 

where Vr{N) = exp{—y^N — lEr/T) is the Boltzmann- 
factor due to the enhanced repulsion from the A — 1 
vortices sharing the columnar defect. Similarly, wn = 
exp(— Z/at-i/T) is the Boltzmann-factor for the bind- 
ing free energy of the A^ — 1 vortices sharing the defect, 
which we approximate by Jn-i ~ Ei or wn « w = 
exp(— Z£'i/r). Considering the generalized exchange be- 
tween these A^ states and noting that exchange of the 
single bound particle is associated with the Boltzmann- 
factor Vexi we arrive at the generalization of 



G{z) 



NG,{z) 



1 - {N - l)VexG,{z) 



The transition point for the exclusion of one particle 
from the exchange is given by the relation Vex,c{N) ~ 
[1 — Vr{N)]/{N — 1). As Vex,c{N) dccrcases for increas- 
ing A^, also states with A^ — 1 or less exchanging par- 
ticles become unstable for Vex < Vex,ciN). This leads 
to the conclusion that particle exchange entirely stops 
at this point, and we thus identify Vex — Vex.c{N) as the 
exchange-delocalization transition point of the A^-particle 
system. For N = 2 our result reduces to Eqs. 1)11(1 and 
(|12|l . whereas we find 



T*\n^/^ (k/IuN) for InN < k 
T*K/\nN for InN > k 



(13) 



for large A^ — B/Bq>. Note that our approach is limited 
to the regime of large vortex spacing a A or A^ ^ 
ajj/X^. The order of the transition is the same as for 
N = 2; i.e., the transition is of infinite order for d ~ 2. 

Phase diagram. In real vortex systems our results 
hold for the limit 5$ « 0. This implies also a vanishing 
vortex density B ^ ii N = B/B^ is fixed. A finite 
vortex density p > corresponds to a system of the fi- 
nite size (X Xj p^l"^ which has no genuine phase transition. 
We thus conclude that the delocalization crossover line of 
Ref. 'H] terminates in the genuine exchange-delocalization 
transition point at i? sa 0, which is given by Eqs. (|12|l or 



(|13ll . For macroscopically large A^ the approximation of 
localized, well-separated exchange and repulsion events 
will break down and our low-density approach will be- 
come invalid whereas the description by a condensate 
of bosonic particles used in Ref. works increasingly 
well in the high-density regime. The exact form of the 
crossover between both descriptions is an open question. 
In both descriptions the delocalization line Bde{T) drops 
exponentially with temperature [see ((13|l ] such that it in- 
tersects with both the pristine melting line i?,„(T) and 
the Bose glass melting line Bbg{T); see Fig. |21 Beyond 
the delocalization line, vortex line exchange at the defects 
sets in, which leads to line wandering and a liquidlike be- 
havior even in the presence of columnar defects, which 
become irrelevant. Therefore, the relevant melting line is 
the pristine melting line for T > Tde{B) in the delocal- 
ized phase, whereas it is the Bose glass melting line in the 
localized phase. Therefore, there exists a range of mag- 
netic fields where the vortex lattice melts by undergoing 
the delocalization transition as shown in Fig. El The 
resulting phase diagram is in good agreement with the 
experimental results regarding the melting of "porous" 
vortex matter 

Conclusion. In conclusion, we have shown that the 
competitive localization of particles with mutual short- 
range repulsion by a short-range attractive defect leads to 
the existence of a genuine phase transition, the exchange- 
delocalization transition, which marks the onset of par- 
ticle exchange at the defect. We have investigated the 
exchange-delocalization transition in a system of vor- 
tices with columnar defects or interacting bosons with 
localizing defects at low density. We expect exchange- 
localization transitions to play an important role in vari- 
ous other systems where competitive localization occurs. 
In the introduction, we have already mentioned the quan- 
tum mechanical system of two interacting electrons com- 
peting for an attractive potential that could be real- 
ized by a quantum dot. We also expect competitive lo- 
calization to be relevant for biopolymers, for example it 
applies to the competitive binding of two identical single 
strands of DNA to a single complementary DNA strand 
as it is important for DNA microarray engineering. 
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